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Abstract 

We study the gravitational scattering of massive particles with and 
without spin in the effective theory of gravity at one loop level. Our 
focus is on long distance effects arising from nonanalytic components 
of the scattering amplitude and we show that the spin-independent 
and the spin-dependent long range components exhibit a universal 
form. Both classical and quantum corrections are obtained, and the 
definition of a proper second order potential is discussed. 
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1 Introduction 

The gravitational interaction of two massive particles is described by New- 
ton's law 



which is a good approximation for the motion of nonrelativistic particles at 
large separations. Einstein's theory of general relativity, however, predicts 
important corrections that have been verified experimentally, e.g., in mea- 
surements of the precession of the perihelion of Mercury pQ. This effect 
can be calculated from the famous Einstein-Infeld-Hoffmann Lagrangian [2] 
where the terms beyond Newtonian physics stem from the relativistic 0(v 2 ) 
corrections to the kinetic energy, from relativistic 0(v 2 ) corrections to the 
Newtonian potential and from a new piece of the potential proportional to 
G 2 M 3 /r 2 which may be regarded as an 0(GM/r) correction to the leading 
Newtonian potential. The two small expansion parameters here are v 2 and 
GM/r which for bound states are related due to the virial theorem. We, 
however, will examine scattering processes wherein the classical expansion 
parameters v 2 and GM/r are independent, and our focus will be on the 
components that vanish in the nonrelativistic limit v — *■ 0. 

In the weak field limit gravitational dynamics can be described in terms 
of a quantum field theory which is based on expanding around flat Minkowski 
spacetime and quantizing the metric fluctuation in terms of a massless spin-2 
field — the graviton. The resulting theory is an effective field theory whose 
predictions are trustworthy only for energies much smaller than the Planck 
scale. Nevertheless, one can in this picture calculate well-defined quantum 
predictions as well as classical corrections to Newtonian physics. Donoghue 
has pioneered the use of the effective field theory of gravity to extract the 
leading long distance effects in the nonrelativistic limit [3l H] (see [5] and [6] 
for reviews) which yield corrections of the form 



where Ac and Aq are the coefficients of the classical and quantum correc- 
tions respectively and are evaluated below. Note that while the classical 
expansion parameter GM/r ~ M/Mpi x ipi/r can give measurable effects 
for macroscopic objects (where a large factor multiplies a tiny factor), the 
quantum corrections which scale as Gh/r 2 ~ (£pi/r) 2 are clearly tiny and 
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phenomenologically irrelevant for macroscopic systems. The calculations of 
these leading long distance corrections are performed by evaluating the lead- 
ing nonanalytic components of the scattering amplitude at one loop level, 
where 0(G 2 ) effects first arise, and by defining a potential in terms of the 
Fourier transform of this amplitude. 

Once spin is introduced into the calculation, additional interaction struc- 
tures arise, such as a spin-orbit coupling and a spin-spin coupling. The 
leading effects of order G stemming from the tree level one-graviton exchange 
process give rise to the familiar geodetic precession and to the Lense-Thirring 
effect [JJ. At the two-graviton exchange level, classical 0(GM/r) and quan- 
tum 0(Gh/r 2 ) corrections to these leading spin-dependent interactions arise 
and are calculated as part of our work. 

In earlier work on loop corrections to the form factors of graviton cou- 
plings it was found that the long distance corrections to the spin-independent 
interactions have the same form for scalars, spin-1/2 fermions and for spin-1 
bosons [3 IE]- These NLO form factor interactions constitute one compo- 
nent of the calculation of the full scattering amplitude at the two-graviton 
exchange level, so the question arises as to whether the full scattering am- 
plitude also exhibits such universalities, whereby the form of the corrections 
to the spin-independent Newtonian potential is independent of the spins of 
the scattered particles — in other words, are the coefficients Aq and Aq in 
Eq. ([2]) independent of spin? 

This then is the goal of the present work — to evaluate the gravitational 
scattering of two massive particles having various spins, in order to check 
previous work and to verify the universality hypothesis. In the next section 
then we review our calculational methods and reproduce the long range grav- 
itational scattering amplitude in the case of a pair of spinless particles. In 
the following sections, we generalize these methods to the cases of spin-0 - 
spin-1/2, spin-0 - spin-1 and spin-1/2 - spin-1/2 scattering and demonstrate 
both the expected universality as well as novel spin-dependent interactions. 
We summarize our work and draw general implications in a brief concluding 
section. Appendix |A] encapsulates parts of our calculational methods while 
we refer to our companion paper on long distance effects in electromagnetic 
scattering [9] for many of the details such as the required Feynman integrals, 
the Fourier transformations etc. In Appendix [B] we demonstrate how the 
classical equations of motion in the form of the Einstein-Infeld-Hoffmann 
Lagrangian can be extracted from our results. 
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a b 

Figure 1: Basic kinematics of gravitational scattering. 

2 Spin-Independent Scattering 

We first set the kinematic framework for our study. We consider the gravi- 
tational scattering of two non-identical particles — particle a with mass m a , 
and incoming four- momentum pi and particle b with mass ra&, and incoming 
four- momentum p 3 . After interacting the final four-momentum of particle a 
is p 2 = Pi — q and that of particle b is p 4 = p 3 + q — cf. Fig. HJ Now we need 
to be more specific. 

2.1 Spin-0 — Spin-0 Scattering 

We begin by examining the gravitational scattering of two spinless particles. 
The gravitational coupling of a spin-0 particle is found by expanding the 
minimally coupled scalar field matter Lagrangian 

V=gC m = (~<r<9 M <R0 - \m 2 <P 2 ) (3) 

in terms of the gravitational field h^ v which is a small fluctuation of the 
metric about flat Minkowski space defined as 

9^u = T]ixv + Kh^ u (4) 

with k = v327rG oc 1/Mp. The inclusion of this factor k in the definition of 
the graviton field h„ u gives this field a mass- dimension of unity and thus yields 
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a kinetic term of standard normalization without a dimensionful parameter. 
For matter interactions, this choice is convenient since the order of k keeps 
track of the number of gravitons involved in an interaction. Once the action 
is written in terms of the expansion of the graviton field, all indices are 
understood to be lowered or raised using the Minkowski metric rj^. We also 
require the expansion of the inverse metric and square root of the determinant 
of the metric tensor — 



g = l + ^h+^(h 2 -2h tIu h^) + 0(K 3 



(5) 



Then, expanding in powers of k, we find 

1„ , 1 



?4? 



-d^d^ - -m 2 2 



-d, 



K 

y 



±(h 2 - 2h lu/ h^) (^d^d^ - \m 2 ^ 
+ (h» a K-^hh^dv 



(6) 



where h = r] a/3 h a p represents the trace and the one- and two-graviton vertices 
are identified as 



i 
i 

P2 k 

I 
I 
I 

I 
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—IK 




PiixPiv + PivP2/i ~ V^iPi -P2-m 2 ) 
1I^ M I^, P MP2+PiP2) 



- ijl^h^pa + V p <jIk\,^)p1P2 ~ P^pcj (pi -P2-m 



(7) 
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where we have defined 

-Pa/3,7,5 = 7^(VayVl3S + VaSVll-y ~ VapV-ys)- ( 8 ) 

The purely gravitational dynamics are derived from the Einstein-Hilbert 
action 

where higher dimensional operators such as R 2 , which are needed when per- 
forming renormalization, are not relevant here [31 H]- Since general relativity 
is invariant under local coordinate transformations, it is a gauge theory and 
one must deal with the subtleties that arise in the quantization of gauge 
theories — we have to perform gauge fixing. The procedure works analogously 
to gauge fixing in Yang-Mills theories. Furthermore, we will make use of the 
background field method, which provides a powerful organizational scheme 
for the quantization of the effective field theory of gravity: Keeping the back- 
ground metric general instead of restricting to the flat Minkowski metric, we 
use the expansion 

g»v = g^v + Khjju/. (10) 

where g^ u is the classical background metric (or field) and h^ v is the quan- 
tum field. A gauge fixing condition is only imposed on the quantum field h^, 
leaving the general covariance of the background unaffected. This procedure 
has the advantage of ensuring that the resulting theory can be renormal- 
ized, since the loop expansion then has the same symmetry properties as the 
action. 

Moreover, the background field method greatly simplifies calculations in- 
volving graviton loops. Gravitons running in loops must be derived from 
from an expansion involving the quantum field whereas gravitons that 
are not within a loop may be derived from expanding the background field 

9fiv = V^u + kH^ (11) 

where denotes an "external" graviton, i.e., a graviton that is not inside 
a loop. At the one loop level, at most two gravitons involved in a vertex are 
propagating within a loop, which allows us to use a gauge fixing condition 
linear in the quantum field h^ u and greatly simplifies the derivation of both 
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the triple graviton vertex and the vertex that couples one graviton to the 
ghost fields. 

In order to fix the gauge, we will use the harmonic gauge — 

- ^Dph = (12) 

where denotes the covariant derivative on the background metric. This 
condition leads to an additional gauge fixing piece of the action 

S GF = J d 4 x^ [p^h^ - l -D^h)j (p p WO - ^D'h) (13) 

as well as the ghost action 

S G host = J d*x^j ff (ppD v - R^) rf (14) 

where rj^ is the ghost field that annihilates a ghost particle while creates 
a ghost particleJll 

Now we are in the position to derive the Feynman rules for the effective 
field theory of gravity. The complete quantum gravitational action consists 
of three components 

SGrav = SqR + SqF + Schost, (15) 

and we will illustrate the use of the background field method during the 
derivations of the Feynman rules. When deriving the graviton propagator 
we expand the action to second order in the quantum fields and to zeroth 
order in the external gravitons H^ u and the ghosts, yielding 

DapM = (16) 

We obtain the ghost propagator expanding the action to second order in the 
ghost fields and to zeroth order in the gravitons, whereby 

DM = (17) 

Our triple graviton vertex is obtained by expanding to second order in the 
quantum fields h^ u , to first order in the background gravitons H^ u , and to 
zeroth order in the ghost fields and reads 

1 Note that the ghost fields anticommute since they obey Fermi-Dirac statistics and we 
have to include a factor of (—1) for each closed ghost loop. 
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T. 



fj,u 



q/3,75 



(k,q) 



af3,-y8 



Wk p + (k + qY{k + q) u + q»q v - -rj^q 2 



rcrv, 



T&u, TXrt. 
1 a/3 1 7 5 



rAcr, r/iZ/, 1 rAcr. TUf, 
1 af3l yS + 1 ryS- 1 a(3 

-q 2 (j] a/3 I^' l5 + rjjsl^'ap) ~ rf v q X q a {r] c ,fil 1 &M + %<s^a/3,Acr] 
+ [2g A (r% 5 / a/3iA(7 P + r^ 1& i a p M k v 

-I aV 'a P IySM(k + qY - I ati 'apI 7 5,Xa(k + qf) 



+v" u q x qA^,x P i pcT ^s + i^x^' 



a/3; 



+ 



(k 2 + (k + q) 2 )[r^ a pl l5 / + I^a/Sl^ ~ ^PaPnS 



{{k + qf Va ^\ 6 + k\ 5 I^ 



a/3) 



(18) 



where the graviton with Lorentz indices [iv represents a background gravi- 
ton, and therefore, it is not to be used within any loop@ The final missing 
Feynman rule from the purely gravitational action is the vertex coupling one 



2 For each closed graviton bubble loop, we also have to include a symmetry factor of 
1/2! in the amplitude. 
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graviton to two ghost fields which we derive by expanding to second order in 
the ghost fields and to first order in H^ v — 



k + q 



a 



IK 



= - (k 2 + (k + q) 2 + q 2 )I a ^ + 2r ]Q pk\k + q yP X(T , 



fj,i/ 



+ 2q a k X Ip X ^ - 2qp(k + q) X I aX ^ + q a q m 



(19) 



where the graviton with Lorentz indices \iv is again a background graviton 
and not to be used as part of any loop. 

We begin our calculation with the simplest tree level single-graviton ex- 
change which leads to the amplitude^! 



°M {1 \q) 



\/2E\2E22E^2E^ 



rll a) (p 2 , Pi) — ^— r ( a/3 (p 4 , Pa) 



q 



-8ttG 



\f2E\2E22E^2E~^ 



m 2 a -m 2 + \q 2 ) 2 - 2m 2 a m 2 b - \q 



2™2 _ 14 
4 



(20) 



A convenient way to define the nonrelativistic potential is as the Fourier 
transform of the nonrelativistic center of mass scattering amplitude. We 
utilize a symmetric center of mass framed with incoming momenta p\ = 
p + q/2 and ps = —p\ = —p — q/2 and with outgoing momenta f>2 = p — q/2 



3 Our normalization of the amplitude is a nonrelativistic one such that after applying 
all Feynman rules we divide the amplitude by a factor of s/2E\ 2Ei2E%2E/^. 

These symmetric momentum labels of the center of mass frame are chosen so that 
the leading order coordinate space potential is real in the calculation of spin-0 - spin-1 
scattering presented below. 
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and p 4 : 



q 2 -. 



p 2 <C m 2 



p + q/2. Conservation of energy then requires p ■ q = so that 
1, 2, 3, 4 and q 2 = —q 2 . In the nonrelativistic limit— 
the lowest order amplitude reads 



p" + q 2 /4: for i 



°M W (q) 



<iixGm a m b 



q 2 



1 + 



P' 



1 + 



+ 



m a m b 



+ Gtt 



3(m 2 a + m 2 ) 
2m a m b 



+ 



p' 



m a m b 



3(m a + m b ) 
2m a m b 



' g _ 5(m 2 a + m 2 ) 2 \ + 
Am 2 a m 2 J 



+ 



(21) 



yielding the potential 

ot/(i)/^ f d3 Q 



(2?r) 3 
Gm a m b 
r 

+ G7i5 3 (r) 



°M {1 \q) e 
p 2 



iqr 



1 + 

m a m b \ 
' 3(m 2 a + m 2 b ) 
2m a m b 



1 



3(m a + m b ) : 



+ . 



+ 



P 



2m a m b 
2 ( 5(ml + ml) 2 " 



m a m b 



4m 2 m? 



+ ... 



(22) 



The leading component of Eq. (122]) is recognized as the usual Newtonian 
potential (accompanied by a small kinematic correction) while the second 
piece is a short range modification. 

Our purpose in this paper is to study the long distance corrections to this 
lowest order potential which arise from the two-graviton exchange diagrams 
shown in Fig. [2j This problem has been previously studied by Iwasaki using 
noncovariant perturbation theory [10], and by Khriplovich and Kirilin [TH 
[12] and by Bjerrum-Bohr, Donoghue, and Holstein [13] using conventional 
Feynman diagrams. Our approach will be similar to that used in [TT], [T2] 
and [13J. That is, using the ideas of effective field theory, we evaluate these 
Feynman diagrams by keeping only the leading nonanalytic structure in q 2 , 
since it is these pieces that lead to the long range corrections to the potential 
while components analytic in q 2 only yield short range contributions, i.e., 
delta functions or derivatives of delta functions. The leading nonanalytic 
behavior is of two basic forms 



i) terms in 1/ \/—q 2 which are ^-independent and therefore classical 
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^=^f^^<> -f- <p^«^^ lp*J^f*m 



(e) (/) (g) 

Figure 2: One loop diagrams of gravitational scattering. 



ii) terms in log — q 2 which are ^.-dependent and therefore quantum me- 
chanical. 

The former terms, when Fourier transformed lead to corrections to the non- 
relativistic potential of the form V c i assica i(r) ~ 1/r 2 while the latter lead to 
y quantum if) ~ h/mr 3 corrections. For typical masses and separations the 
quantum mechanical forms are themselves numerically insignificant. How- 
ever, they are intriguing in that their origin appears to be associated with 
zitterbewegung. That is, classically we can define the potential by measuring 
the energy when two objects are separated by distance r. However, in the 
quantum mechanical case the distance between two objects is uncertain by 
an amount of order the Compton wavelength due to zero point motion — 
5r ~ h/m. This leads to the replacement 

1 1 1 % 



3 



r 2 (r ± 5r) 2 r 2 mr 

which is the form found in our calculations. 

The diagrams to be evaluated are shown in Fig. [2] where the "blobs" 
shown in diagrams (f) and (g) are explained in Fig. [3j The calculational 
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-\- «*s» f*> 

Figure 3: Loop corrections subsumed in vertex and in vacuum polarization 
functions. 



details including the Feynman integrals are described in an appendix of our 
companion paper on electromagnetic scattering [9] . Here we present only the 
results. Dennin£n 



S 



7[- 



and 



log -q 2 



we have, from diagrams [2](a)-(g) respectively 



G 2 m a m b [-44L] 
G 2 m a m fe [28L + 8m a S] 
G 2 m a m b [28L + 8m fe S] 

G 2 m a m b 



4m a m b 8(m 2 a + m 2 ' 



iAnG 2 m 2 mi — 
q 2 



lm a m b 



G 2 m a m b 



s - s 
4m a m 6 



m a m b 



\(m 2 a + m 2 b ) 
m a m b 



8}L + 4:(m a + m b )S 



70 
~3~ 



L -4(m a + m b )S' 



5 Note that our definition of S differs slightly from previous publications [3J [4j [T3l [T] in 
that it does not include a factor of mass in the numerator. 
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G 2 m a m b [14L - 2(m a + m b )S] 
G 2 m a m b 



— -L 

15 



(23) 



where s = (pi+ps) is the square of the center of mass energy and so = (m a + 
irib) 2 is its threshold value. The contributions from the form factor diagrams 
in Fig. E](f) have been calculated both directly and using the results from 
[3 E] , and the contribution from vacuum polarization diagrams in Fig. g) 
have been obtained both by direct evaluation and by using previous results 
of 't Hooft and Veltman [TJ] for the divergences of the graviton self-energy 
that allow us to infer its nonanalytic contribution^]. Here we should note 
that any massless particle would contribute to the nonanalytic terms of the 
graviton vacuum energy since gravitons couple to anything. We choose here 
not to include other massless particles besides the graviton — in particular we 
do not include the contribution of the photon in the quantum piece of our 
potential. 

Summing, we find the total result 



M { $ t {q) = G 2 m a m b 



41 ' 

6(m a + m h )S - —L 
5 



iAjcG 2 m 2 ml — 1 



im a m b 



(24) 



s - s 

and we observe that, in addition to the expected terms involving L and S, 
there arises a piece of the second order amplitude which is imaginary. The 
origin of this imaginary piece is, of course, from the second Born approxi- 
mation to the Newtonian potential, and reminds us that in order to define 
a proper correction to the first order Newtonian potential we must subtract 
off such terms. For this purpose we will work in the nonrelativistic limit and 
the center of mass frame — p\ + p% = as defined above. We have then 



s - s = 2 J ml + piJm 2 + p{ + 2p{ - 2m a m b 



and 



im a m b 



V s - so Po 

where m r = m a m b / {m a + m b ) is the reduced mass andpo = \Pi\, i 
The transition amplitude Eq. (j24"l) then assumes the form 



(25) 



(26) 



1,2,3,4. 



°A^S(g)~G 2 m a m 6 



41 

6(m a + m b )S - —L 
5 



- l ^G 2 m 2 a m 2 b ^. (27) 

<r Po 



6 In a massless theory the divergences of dimensional regularization are accompanied 
by logarithms of the momentum transfer. 
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For the iteration we shall use the simple potential 



Gm a irib 



(28) 



which reproduces the long distance behavior of the lowest order amplitude 
for spin-0 - spin-0 gravitational scattering — Eq. fl22|) — in the nonrelativistic 
limit. We will employ the momentum space representation 



where we identify pi 
°Ampg ) (g) 



Pf 



Pi 



q 2 



(Pi - Pf) 2 

Pi and pf = P2, and the second Born term is then 

dH 



(29) 



Pf 



V G 



0t/(1) 



P, 



(2tt) 3 
dH 



E( Po ) - E(£) 
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where 



(2tt) 3 
G (0) (£) = 



Pf) g (0) $ °y G l) - n (so) 



(31) 



+ ie 



2m r 2m r 

is the nonrelativistic propagator. Note that in Eq. (I30I) we take both the 
leading order potential as well as the total energies E(p ) and E(£) in the 
nonrelativistic limit. The remaining integration can be performed exactly, 
yielding^ 



d 3 £ —A^Gm a mb 



-47tGm a mb 



(2tt) 3 \£ 



p*\ 2 + x2 £t-£;.+^ I* 



+ A 2 



H 



rj r\ r\ T 1 TYXiy 

-iAnG m a m b — — 
T Po 



(32) 



which precisely reproduces the imaginary component of °Mtot(<l)) as ex_ 
pected. In order to produce a properly defined second order potential Vq^ (r) 
we must then subtract this second order Born term from the second order 



r Note that the iteration integrals are listed in Appendix [A] 
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transition amplitude, yielding the result 



' f '''< -" r, ' ro Xg(9)-°Ampg»( g ) 



G m a m b 



41 ' 

-6S(m a + m b ) + —L 



(2tt) 3 

3G 2 m a m b (m a + m b ) A\G 2 m a m b h 



5 



r 2 107rr 3 



(33) 



The quantum mechanical — ~ h/r 3 — component of the second order po- 
tential given in Eq. ( 1331) agrees with that previously given by Bjerrum-Bohr, 
Donoghue, and Holstein [13] and by Kirilin and Khriplovich [12]. However, 
the classical — ~ 1/r 2 — contribution quoted by Iwasaki 



o T /(2)/^ _ G 2 m a m b {m a + m b ) 



W(r) = " 2r2 " (34) 

differs from that quoted above in Eq. (|33|) and by Bjerrum-Bohr et al. in [T3] . 
The resolution of this issue was given by Sucher, who pointed out that the 
classical term depends upon the precise definition of the first order potential 
used in the iteration [15]. Moreover, it depends and on whether one uses 
relativistic forms of the leading order potentials and the propagator G^ (£) in 
the iteration. In modern terms, the potential depends on how one performs 
the matching — e.g., Iwasaki [10] performs an off-shell matching while we 
match on-shell. In our companion paper on electromagnetic scattering [9] we 
have provided a more detailed discussion of these ambiguities^. Use of the 
simple lowest order form Eq. (I2TJ1) within a nonrelativistic iteration yields our 
result for the iteration amplitude given in Eq. (132]) and is sufficient to remove 
the offending imaginary piece of the scattering amplitude. In Appendix [B] 
we derive an alternative form of the 0(G 2 ) classical potential which results 
from an iteration that includes the leading relativistic corrections and which 
reproduces the classical equations of motion. 

Therefore, a unique definition of the potential does not exist. Of course, 
the ambiguities in the form of the second order classical potential should not 



8 Besides the dependence on the forms used in the iteration, the classical piece also 
depends on the coordinates used. The quantum piece however depends neither on the 
choice of coordinates [13] nor on the iteration forms [9] . 
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be a concern, since the potential is not an observable. What is an observable 
is the on-shell transition amplitude, which is uniquely defined in each case 
as 

+ °Am ft (?1 (35) 

where the index % denotes differing possible definitions of the potentials and 
the iteration. Thus we regard our potential as a nice way to display our re- 
sulting scattering amplitudes in coordinate space, but we emphasize that our 
main results are the long distance components of the scattering amplitude. 



°Mtot{q) = ~ 



d 3 re^ 



Spin-Dependent Scattering: Spin-Orbit In- 
teraction 



3.1 Spin-0 - Spin-1/2 

Having determined the form of the potential for the spinless scattering case 
we move on to the case of scattering of particles carrying spin. We begin 
with the scattering of a spinless particle a from a spin-1/2 particle b. 

For the case of spin-1/2 we require some additional formalism in order 
to extract the gravitational couplings, which is necessary because the Dirac 
algebra {7 a ,7 b } = 2i] ab is defined with respect to the Minkoswki flat space 
metric. In this case the Dirac matter Lagrangian coupled to gravity reads 



1> 



(36) 



and involves the vierbein e M a which links global coordinates with those in 
a locally flat space. The vierbein is in some sense the "square root" of the 
metric tensor g^ v and satisfies the relations 



e a e 



Vab = 



b UU 

v y 



Vab 



e\e v brf 

^ 1 a 6 b 9 pi/ 



= 9 
= V 



ab 



The covariant derivative is 

D, 



df: R + 



,Vac 0- 



cb 



(37) 



(38) 



with a cb 



7,7 



in such a way that 



2 1 4 ~M b' 

and the partial derivative d^ R acts only on spinors and 



(39) 
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Putting everything together, we find then 



m 



ip. (40) 



The spin connection uj fl a b rj ac can be derived in terms of vierbeins by requir- 
ing D^e v a = and by antisymmetrization in /x <-> v in order to get rid of 
Christoffel symbol^. The result is: 



Vaf 



b^c) (41) 



In order to derive the Feynman rules we expand the ingredients in Eq. ( 1401 ) 
that contain graviton couplings, that is we need e^ a and Up a b r] ac expanded 
up to 0{k 2 ) 



K k/^ 

+ — h p h dph cp - — h p b d p hp C + — h p b d c hp P 
2 o 4 4 



(42) 



After these expansions are employed, we no longer need to distinguish be- 
tween Latin Lorentz indeces and Greek covariant indices and can use the 
Minkowski metric to lower and raise all indices. 

The matter Lagrangian then has the expansion — (note here that our con- 
ventions are 75 = — Z7°7 1 7 2 7 3 and e 0123 = +1) 



-7(1 



4 



— m )ip 



(i) 

m 



9 For our purposes we shall use only the symmetric component of the vierbein matrices, 
since these are physical and can be connected to the metric tensor, while their antisym- 
metric components are associated with freedom of homogeneous transformations of the 
local Lorentz frames and do not contribute to nonanalyticity [16] . 
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"iff 
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2 ( 



(43) 



and the corresponding one- and two-graviton vertices are found to be 
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Pi JL 
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Vpu(lp(Pl +P2)a + l*(Pl +P2)p) 
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(44) 

The tree level transition amplitude from one-graviton exchange is then 
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Defining the spin vector as 



we find the identity 
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whereupon the nonanalytic part of the transition amplitude in the thereshold 
limit s — > s — ( m a + ^f>) 2 can be written in the form 
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In order to define the potential we require the nonrelativistic amplitude in 
the symmetric center of mass frame (pi = — p 3 = p + q/2) where 
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and the lowest order potential becomes 
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r 1 r 6 2m b 
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Gm a m b 



(53) 



where L = f x p is the angular momentum — the modification of the leading 
spin-independent potential has a spin-orbit character. 

In order to determine the second order potential, we must evaluate the 
one loop diagrams in Fig. [2j A subtlety that arises in the calculation involv- 
ing spin is that two independent kinematic variables arise: the momentum 
transfer q 2 and s — Sq, which is to leading order proportional to p\ (where 
Pq = p 2 , i = 1,2,3,4) in the center of mass frame. We find that our re- 
sults differ if we perform an expansion first in s — sq and then in q 2 or vice 
versa. This ordering issue occurs only for the box diagram, diagram (d) of 
Fig. [21 where it stems from the reduction of vector and tensor box integrals. 
Their reduction in terms of scalar integrals involves the inversion of a ma- 
trix whose Gram determinant vanishes in the nonrelativistic threshold limit 
q 2 , s — So — > 0. More precisely, the denominators or the vector and tensor box 
integrals (see Appendix A in [9]) involve a factor of (4po — <f 2 ) when expanded 
in the nonrelativistic limit. Since q 2 = 4p§ sin 2 | with 9 the scattering angle, 
we notice that 4p\ > q 2 unless we consider backward scattering where 9 = n 
and where the scattering amplitude diverges. And since p^ originates from 
the relativistic structure s — Sq, it is clear that we must first expand our 
vector and tensor box integrals in q 2 and then in s — sq. With this procedure, 
the contributions of the loop diagrams in Fig. [2] to the spin-0 - spin-1/2 
scattering amplitude read 



G 2 m a m b u(p 4 )u(p 3 )(-10L) H u(p A ) ^i«(p 3 )(-24L) 
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(54) 



Again, we have calculated the form factor diagrams in Fig. [2](f ) both directly 
and via the results from [8] . Summing, we find 
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Using the identity Eq. (j4"7j) and 



Pi ■ (P3 + Pa) = 2m a m b + s - s + — 



Eq. (1551) becomes 
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(56) 



Finally, working in the center of mass frame and taking the nonrelativistic 
limit, we find 
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-S b -pxq (57) 



We note from Eq. (1571) that the scattering amplitude consists of two pieces — a 
spin-independent component proportional to XjX\ whose functional form 



/ 41 \ 
G 2 m a m b f 6(m a + m b )S —L J — iAnG 2 m 2 a m 2 



L m r 
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(58) 



is identical to that of spinless scattering — together with a spin-orbit compo- 
nent proportional to 

— S b -p x q 
m b 
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whose functional form is 
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(59) 



We note in Eq. (1591) the presence in the spin-orbit potential of an imaginary 
final state rescattering term proportional to i/po, similar to that found in 
the case of spin-independent scattering, together with a completely new type 
of kinematic form, proportional to 1/Pq which diverges at threshold. The 
presence of either term would prevent us from writing down a well defined 
second order potential. 

The solution to this problem is, as before, to properly subtract the iterated 
first order potential — 
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where we now use the one-graviton exchange potential ^Vq\t) given in Eq. 
( I5"3"|) . Splitting this lowest order potential into spin-independent and spin- 
dependent components — 
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(62) 



we find that the iterated amplitude splits also into spin-independent and 
spin-dependent pieces. The leading spin-independent amplitude is 
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4:7rGm a mb, and the leading spin-dependent term is 
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(In principle we would also have to iterate the leading order spin-orbit piece 
twice. However this procedure yields only terms higher order in p 2 ,.) We 
observe that when the amplitudes Eqs. (l64"j) and fl63l are subtracted from 
the full one loop scattering amplitude Eq. (j57p both the terms involving 
1/Pq and those proportional to i/po disappear leaving behind a well-defined 
second order potential 

d 3 q 



-iqr 



MS(g)-fAmp^(g) 



(2), 



d 3 q 
(2tt)2 



-iq-r 



41 



G m a m 6 I -6(m a + m^S 1 + — L J \]x 



G - 



12m 3 + 45m 2 m& + 56m a m 2 + 24m 3 , 
2(m a + m b ) 



S 



23 



+ 



, 87m a + 128m b \ i g _ 

+ T7\ L — S b -pxq 

10 / m h 



3G 2 m a mb(m a + mj) 41G 2 m a mbh 



107rr 3 



b] b 

XfXi 



1 Pt -. A 

+ — S b -pxV 



G 2 (12m 3 + 45m 2 ! m fe + 56m a m^ + 2Am\ 

4(m a + mb)r 2 
G 2 (87m a + I28m b )h 



20nr 3 

3G 2 m a mb(m a + m&) 41G 2 m a mbh 



r 2 107rr 3 
G 2 {12m? a + Ahm 2 a m h + 5Qm a ml + 24ml 

2mb(m a + rrib)r 4 
3G 2 (87m a + 128m b )h 



b] b 
Xf Xi 



20nmbr 5 



L-Sb 



(65) 



We observe then that the second order potential for long range gravitational 
scattering of a spinless and spin-1/2 particle consists of two components: one 
which is independent of the spin of particle b and is identical to the potential 
found for the case of spinless scattering, accompanied by a spin-orbit inter- 
action involving a new form for its classical and quantum potentials. It is 
tempting to speculate that the form of this new spin-orbit potential is also 
universal. In order to check this hypothesis we consider the case of spin-0 - 
spin-1 scattering. 



3.2 Spin-0 - Spin-1 



The dynamics of a neutral spin-1 field M having mass m is described by the 
Proca Lagrangian which, when coupled to gravity via minimal substitution, 
takes the form 
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where 



U, 
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0, 



<%0 M (67) 

is the spin-1 field tensor. The last equality in Eq. fl67|) follows from the 
symmetry of the connection coefficients r" = r™ . Expanded in terms of 
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the graviton field, the matter Lagrangian then has the form 
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and the one- and two-graviton vertices are 
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(69) 

If we take the incoming spin-1 particle to have polarization vector e b sat- 
isfying • P3 = and the outgoing particle to have polarization e b satisfying 
e b f -P4 = 0, then the one-graviton exchange amplitude can then be written as 
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where the first expression is exact while the second expression contains only 
the nonanalytic part in the threshold limit s — > s = ( m a + rnb) 2 ■ Now we 
rewrite this expression using the identity 



e fp e i Q e ip, e f Q — 



4ml 



(3 -y riS 

tp^sPsq o b 



(P3+P±)p m a( p _ 
2ml f q 1 q 



(71) 



26 



where we have defined the spin vector 



The leading one-graviton exchange amplitude can then be written as 
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The spin-independent and spin-orbit terms here are identical in form to those 
found in the spin-0 - spin- 1/2 case but now are accompanied by new terms 
which are quadrupole in nature, as can be seen from the identity 
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The corresponding lowest order potential is then 
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where we have defined 

w : T b : s = w c T h cd s d 

and which agrees precisely with its spin- 1/2 analog in Eq. (|53|) up to qua- 
drupole and tensor corrections. 

The calculation of the one loop corrections proceeds as before, but with 
increased complexity due to the unit spin. Evaluating the diagrams in Fig. 
HI we find then 
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Combining, we find the full one loop amplitude 
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which, using the identity Eq. ( 17T1) . becomes 
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Notice here that without the e b f * ■ qe b ■ q terms, Eq. 



has an identical 



structure to that of the case of spin-0 - spin-1/2 scattering — Eq. (1551) - 
provided we substitute u{pi)u{pz) - 



-> e f e i- 

Finally, taking the nonrelativistic limit we find 
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As found in the earlier calculations, there exist terms involving both i/po 
and 1/po which prevent the defining of a simple second order potential. The 
solution now is well known — subtraction of the iterated first order potential. 
Since the form of the spin-independent — e* B ■ — and spin-orbit — S b ■ p x q- — 
terms is identical to that found for the case of spin-1/2, it is clear that the 
subtraction goes through as before and that the corresponding pieces of the 
second order potential have the same form as found for spin-1/2. In addition, 
there are now two new pieces of the amplitude, the quadrupole structure 
q :T b : q which multiplies terms involving both i/p and 1/pjj and the tensor 
structure p : T b : p multiplying only i/po- In order to remove these we must 
iterate the full first order potential including these quadrupole and tensor 
components. However, we find that our simple nonrelativistic iteration fails 
to remove them! We suspect the reason to be the presence of the tensor 
structure p : T b : p in the lowest order potential which is in some sense a 
relativistic correction but which when iterated yields also quadrupole pieces 
q : T b : q. A fully relativistic iteration should thus be performed which is 
under study but is beyond the scope of this paper. It would be interesting 
to investigate if the requirement of canceling all i/po and 1/pjj forms in the 
quadrupole and tensor pieces could clarify the ambiguity in the iteration of 
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the leading order potential as discussed in [T5J E] • 

Since we did not perform the proper iteration of the quadrupole and 
tensor pieces we include only the spin independent and spin-orbit pieces in 
the resulting second order potential 
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where 1 V^ 2 ' ) (r*) denotes the tensor pieces not explicitly shown. Comparison 
with the corresponding form of 2V^(f^ given in Eq. (165"]) confirms the 
universality which we have suggested — the spin-independent and spin-orbit 
terms have identical forms. The next task is to see whether this universality 
applies when both scattered particles carry spin. For this purpose we consider 
the case of spin-1/2 - spin-1/2 scattering. 



4 Spin-Dependent Scattering: Spin-Spin In- 
teraction 

4.1 Spin-1/2 - Spin-1/2 

In the case of scattering of a pair of spin-1/2 particles, the basic vertices have 
already been developed in the spin-0 - spin-1/2 scattering section, so we can 
proceed directly to our calculation. The one-graviton exchange amplitude at 
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tree level reads 
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for system a, Eq. (186]) can be written as 
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In the symmetric center of mass frame we take the nonrelativistic limit yield- 
ing 
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whereby the lowest order potential for spin- 1/2 - spin- 1/2 scattering becomes 
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Note that since the piece proportional to S a ■ S b in Eq. (189]) is analytic in q 2 
it yields only a short distance contribution which is omitted in the potential 
Eq. pp. 

In this case when we evaluate the loop diagrams in Fig. [2J we notice that 
part of the spin-spin structure piece contains the form q 2 S a ■ S b multiplying 
the nonanalytic structures L and S. Due to the extra factor of q 2 in this 
form, we must expand all loop integrals to one order higher in q 2 than before 
in order to be consistent. This has been done and does make a difference in 
our results for the spin-spin interaction piece. The results for the individual 
diagrams are then found to be 
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with i = a, b to keep our notation compact. The sum of all diagrams is found 
to be 
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Comparing with our finding for spin-0 - spin-1/2 scattering in Eq. (|56|) . 
we notice the by now familiar universality of the amplitude: The form of 
the component proportional to u{p4)u(j) 3 ) of Eq. (151)1) is found here in the 
component proportional to U a U b , and the form of the structure proportional 
to tafrsPiPsQ 1 'Si °f Eq. (156|) is now found in the component proportional 
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to U a Sb in Eq. (I94p . Moreover, the amplitude is symmetric in a «-> b and 
new gravitational spin-spin interaction corrections arise. The quantum part 
of the spin-spin component has been calculated previously by Kirilin [T7] 
whose result disagrees with our result for the numerical prefactors. For the 
quantum terms in the spin-orbit components however, we fully agree with 
Kirilin's results in [TT] . In the nonrelativistic limit we obtain the expression 
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which obviously has to exhibit the universalities of the spin-independent and 
the spin-orbit pieces. Therefore, we know that the subtraction of the second 
Born iteration successfully removes the unwanted i/po and 1/pq structures 
for the spin-independent and the spin-orbit components. 

The leading spin-spin term of the second Born iteration amplitude is new, 
however, and we compute 
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where we again defined Cq = —AnGm a m b . With this, the full second Born 
iteration amplitude becomes 
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and we observe that when this amplitude is subtracted from the full one loop 
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scattering amplitude Eq. (I9l>j) . all terms involving I /pi and i/p Q disappear 
leaving behind a well-defined second order potential 
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which besides the universal spin-independent and spin-orbit components dis- 
plays a new (and presumably universal) gravitational spin-spin interaction. 



5 Conclusions 

Above we have analyzed the gravitational scattering of two particles hav- 
ing nonzero mass. In lowest order the interaction arises from one-graviton 
exchange and leads at threshold to the well known Newtonian interaction 
V{r) = —Gm a mb/r. Inclusion of two-graviton exchange effects means adding 
the contribution from box, cross-box, triangle, and bubble diagrams, which 
have a rather complex form. The calculation can be simplified, however, 
by using ideas from effective field theory. The point is that if one is inter- 
ested only in the leading long-range behavior of the interaction, then one 
need retain only the leading nonanalytic small momentum-transfer piece of 
the scattering amplitude. Specifically, the terms which one retains are those 
which are nonanalytic and behave as either 1/^—q' 2 or log — q 2 . When Fourier 
transformed, the former leads to classical (/i-independent) terms in the po- 
tential of order G 2 M 3 /r 2 while the latter generates quantum mechanical (h- 
dependent) corrections of order G 2 M 2 h/r 3 . (Of course, there are also shorter 
range nonanalytic contributions than these that are generated by scattering 
terms of order q 2n /y/—q' 2 or g 2n log— q 2 . However, these pieces are higher 
order in momentum transfer and thus lead to shorter distance effects than 
those considered above and are therefore neglected in our discussion.) 

Specific calculations were done for particles with spin — 0, — 1/2, — 1, 
and 1/2 — 1/2 and various universalities were found. In particular, we found 
that in each case there was a spin-independent contribution of the form 



where L = log — q 2 and S = n 2 / \J—q 2 and with S a the spin of particle a 
and Sb the spin of particle b with projections m a and m b on the quantization 
axis. The imaginary component of the amplitude, which would not, when 
Fourier-transformed lead to a real potential, is eliminated when the iterated 
lowest order potential contribution is subtracted, leading to a well defined 
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spin-independent second order potential of universal form 
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whose classical component depends on the way the iteration of the leading 
order potential is performed. This ambiguity shows that the second order 
potential itself is not an observable, but we use it as a nice way to display the 
long distance components of the scattering amplitude in coordinate space. 

If either scattering particle carries spin then there exists an additional 
spin-orbit contribution, whose form is also universal 
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where we have defined 



S a — \S a ,m a f S S a ,m ai j and S b = {S b ,m bf S S b ,m b , 

In this case a well defined second order potential required the subtraction of 
infrared singular terms behaving as both i/p Q and l/p$ which arise from the 
iterated lowest order potential. 

In the calculation of spin-0 - spin-1 scattering we encountered new tensor 
structures including a quadrupole interaction. Unfortunately, the subtraction 
of the i/po and 1/pq tensor pieces in the two-graviton exchange amplitude was 
not successful with our simple nonrelativistic iteration of the leading order 
potential so that we cannot at this time give the form of the quadrupole 
component of the potential. Further work is needed to clarify this issue. The 
corrections to the spin-spin interaction have only been calculated in spin-1/2 
- spin-1/2 scattering where we found their contributions to the scattering 
amplitude and to the potential. Since we verified these forms only for a 
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single spin configuration we have not confirmed its universality which we, 
however, strongly suspect. Of course, for higher spin configurations, there 
also exist quadrupole-quadrupole interactions, spin-quadrupole interactions, 
etc. However, the calculation of such forms becomes increasingly cumbersome 
as the spin increases, and the phenomenological importance becomes smaller. 
Thus we end our calculations here. 

One point of view to interpret the universalities of the long distance com- 
ponents of the scattering amplitudes and the resulting potentials is that if 
we increase the spins of our scattered particles, all we do is to add additional 
multipole moments. The spin-independent component can then be viewed as 
a monopole-monopole interaction, the spin-orbit piece as a dipole-monopole 
interaction etc. As long as we do not change the quantum numbers that 
characterize the lower multipoles, an increase in spin of the scattered parti- 
cles merely adds new interactions that are less important at long distances. 
The same kind of universalities were also found in long distance effects in 
electromagnetic scattering [9] and in the long range components of mixed 
electromagnetic-gravitational scattering [TS] . 
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A Iteration Integrals 

In this appendix we give the integrals 
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which are needed in order to perform the iteration of the lowest order Newton 
potentials. Here we list only the results; for a more detailed derivation, albeit 
with a different prefactor, see [9]. The leading expressions for the iteration 
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integrals read 
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B Classical Equations of Motion 

Above we have argued that the scattering amplitude which is ultimately 
related to observables in quantum field theory is a physical result while the 
potential we have given is not an observable and depends both on the gauge, 
i.e., the choice of coordinates used, and on the way we perform the iteration, 
i.e., on the way we perform the matching. While the classical component 
of our potential is in fact plagued by these ambiguities, the quantum part 
is unique since it is unaffected by how we perform the matching and since 
a quantum field theory calculation in any gauge would result in the same 
result [13]. 

In this appendix we will demonstrate how we can recover the classi- 
cal equations of motion from our scattering amplitudes by setting up the 
Einstein-Infeld-Hoffmann (EIH) Lagrangian [2]. The EIH Lagrangian is it- 
self dependent on the choice of coordinates, but can be expressed in the 
center of mass frame (P = p a = —pb, r = r a — r&) in a general way as [19j 120] 



where the kinetic energy to NLO in the nonrelativistic expansion reads 



T-V 



(104) 



T 



p2 p2 p4 p4 



(105) 



2m a 2rrib 8ml 8mf 



and the potential is 



v = v (1) + v {2) 



(106) 
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with 
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The parameter a parameterizes the choice of coordinates used, where a = 
was the gauge of the original EIH result. The coordinate change 



r 1 — a 
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which implies 



P -> P + a 



G(m a + m b ) 



P - P-r)r 



(109) 



110) 



brings the original EIH Lagrangian into the form above, which is the most 
general result. 

Since we perform our matching on-shell, i.e., we use the on-shell one- 
graviton exchange amplitude to define the leading order O(G) potential, 
terms proportional to P ■ f would never arise. Clearly, our result must be in 
a gauge such that the coefficient of the structure 

[P-ff 

m a m b 

in Eq. 11 10 7ft vanishes. That is the case if and only if the gauge parameter is 

m a m b 



a 



2{m a + m b ) 2 
whereby the EIH potential becomes 
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(113) 
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Comparing the EIH potential in this gauge of Eq. (j!12p with the long 
distance component of the leading order spin-independent potential in Eq. 
(122]) we find full agreement for the relativistic corrections to the 0(G) poten- 
tial. However, comparing the EIH potential in this gauge of Eq. (11131) 
with the classical component of our spin-independent potential in Eq. (1331) 
we see that the two do not agree! The reason for this discrepancy is that 
we elected to use a nonrelativistic iteration when we performed the second 
Born iteration of the leading order potential in Eq. (I32I) . This procedure, 
however, is not self-consistent when we are interested in equations of motion 
at NLO, for which we must account for the leading relativistic corrections in 
the iteration. In particular, we must use expressions for the potential and the 
propagator in Eq. (1301) which include the leading relativistic correctional 
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which yields a second Born iteration amplitude 
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Subtracting this iterated amplitude which includes all corrections to NLO 
from the scattering amplitude °Aitot(o) °f Eq. (I2~T|) we find then the second 
order potential 



0t/(2) 
V i 



NLOV) 



d 3 q 



-iq-r 
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10 The subscript NLO in this sections refers to the iteration being performed at NLO in 
the relativistic expansion. 
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and observe that now the classical component agrees with the 0(G 2 ) EIH 
potential of Eq. (11 13|) . 

Thus we have shown that if we consistently take into account the v 2 
and GM/r corrections beyond Newtonian physics we reproduce the EIH 
Lagrangian in a certain gauge. From the resulting EIH Lagrangian we could 
evaluate observables such as the precession of the perihelion of Mercury which 
must clearly be independent of the gauge used. The inclusion of the v 2 
corrections is required since the equations of motion can be used to describe 
bound states where v 2 ~ GM/r by the virial theorem. 

However, our methods are clearly clumsy for the calculation of classical 
observables. Recently, Goldberger and Rothstein have developed an effective 
field theory of gravity which is optimized for calculating classical observables 
of bound states called NRGR [2D [23 123 121 123 126]- Here the external 
particles are static sources so that no loops are to be calculated in their theory 
when calculating classical observables since the only propagating particles 
present are gravitons which are massless and thus the loop expansion in 
NRGR corresponds to an expansion in h. In the NRGR framework the spin- 
dependent classical equations of motion were calculated recently to NLO by 
Porto and Rothstein [23 [23 123 EH] so that we will not continue here to 
evaluate the corresponding spin-dependent classical potentials consistently 
taking into account all relativistic 0(v 2 ) effects in the iteration. 
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